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ABSTRACT 

This  paper  develops  a  new,  more  accurate  temporal  model  of  phase  screen  generation.  The  long  standing  Fourier 
transform  (FT)  based  method  assumes  the  frozen  flow  hypothesis  holds,  where  large  phase  screens  are  generated 
and  then  shifted.  The  result  is  a  statistically  correct  screen.  Realistically  however,  the  phase  changes  with  time, 
especially  when  the  wind  velocity  is  small  or  non-existent.  The  temporal  evolution  method  proposed  in  this  paper  is 
based  on  expanding  a  random  walk  algorithm  to  the  Zernike  polynomial  method  that  will  include  not  only  the 
perceived  shifting  of  the  Kolmogorov  FT  method,  but  also  a  "boiling  effect"  that  changes  the  phase  as  it  shifts.  This 
new  method  of  phase  screen  generation  will  be  validated  through  a  simulated  experiment  which  measures  the 
correlations  of  tilt  as  a  function  of  time  and  compares  that  to  a  predicted  tilt  correlation  function  derived  using  both 
the  frozen  flow  hypothesis  and  the  new  model.  This  validation  experiment  will  show  that  the  frozen  flow  model 
alone  fails  to  accurately  predict  the  temporal  correlation  of  optical  tilt. 

1.  INTRODUCTION 

For  conventional  imaging  systems.  Geosynchronous  Earth  Orbit  (GEO)  space  objects  cannot  be  resolved  due  to 
their  40  Mm  distance.  The  need  to  obtain  high  resolution  images  of  GEO  objects  is  quite  real  and  to  accomplish  this 
task,  the  Air  Eorce  Research  Laboratories  (AERL)  is  currently  investigating  the  suitability  of  inverse  synthetic 
aperture  laser  radar  (ISAL)  at  the  Air  Eorce  Maui  Optical  and  Supercomputing  Site  (AMOS).  A  critical  component 
to  determining  its  suitability  is  to  accurately  model  the  atmospheric  impacts  on  LADAR  pulses.  The  expected  result 
is  that  atmospheric  impacts  are  a  strong  function  of  the  illumination  wavelength,  by  the  conventional  knowledge  that 
shorter  wavelengths  are  more  affected  than  longer  ones.  Conventional  thinking  is  that  while  the  atmosphere  churns, 
wind  is  the  predominant  effect  that  simplifies  all  modeling  and  simulation  into  the  frozen  flow  hypothesis.  The 
concern  is  that  frozen  flow  hypothesis  based  phase  screen  generation  techniques  fail  to  accurately  predict  the 
temporal  correlation  of  optical  phase  in  low  wind  velocity  situations,  thus  the  intent  of  this  work  is  to  develop  a 
turbulent  flow  model.  This  "boiling  effect"  is  the  churning  and  fluctuation  in  realistic  turbulence.  We  derive  the 
temporal  evolution  aspect  from  [4]  who  applies  a  random  walk  to  the  first  two  Zernike  polynomials,  accurately 
simulating  atmospheric  tilt.  Wave  front  phase  aberrations  introduced  by  the  atmosphere  are  dominated  by  tilt, 
causing  the  most  degradation  of  performance  in  direct-detection  LADAR  systems.  This  project  expands  upon  [4]'s 
tilt  algorithm;  and  is  capable  of  including  as  many  N  number  of  Zernike  polynomials  as  desired  for  low  wind 
velocity  environments. 


2.  BACKGROUND 

We  simulate  atmospheric  turbulence  phase  screens  with  MatLab  through  two  common  modeling  approaches: 
modified  von-Karman,  and  Zernike  polynomial-based  phase  screen  generation.  We  initiate  this  investigation  by 
introducing  and  defining  the  Eourier  transform  (ET)  to  which  we  compare  our  new  Zernike-based  method  of 
temporal  phase  screen  generation. 

2.1  Fourier-based  Phase  Screen  Generation 

The  ET  methods  are  most  common  since  very  large  phase  screens  can  be  generated  quickly.  Most  ET  methods  are 
similar  in  that  they  are  based  on  their  respectively  named  power  spectral  density  (PSD).  One  main  downside  is  to 
this  approach  is  that  lower-order  aberrations  such  as  tilt  are  often  under-represented  [2].  These  lower-order 
aberrations  make  up  a  majority  of  the  atmospheric  energy  spectrum  and  as  such,  must  be  included  to  produce 
realistic  models.  One  of  the  more  widely  used  ET  based  methods  is  based  on  the  von-Karman  spectrum. 
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cI>(k)  =  0.023(ro)“^/^e  (k^  + 

The  modified  von-Karman  allows  an  estimation  of  the  inner  (Ig)  and  outer  (Lg)  scales  of  turbulent  eddies,  leading  to 
K„  =  5.92/lg,  and  Kg  =  2nlLg,  in  which  k  is  spatial  frequency  and  is  scaled  by  a  factor  of  Fried's  constant,  rg.  For  this 
investigation  we  approximate  ideal  limits,  Lg  approaches  infinity  and  Ig  approaches  zero.  To  better  reflect  the 
AMOS  observatory,  we  select  D  to  be  1.6  m  and  rg  to  be  14  cm.  With  these  parameters  we  simplify  Eq.  (1)  to  Eq. 
(2). 


<P(k)  =  (2) 

As  ET-based  methods  are  not  the  focus  of  this  paper  we  limit  our  discussion  to  this  point.  Eor  comparative  analysis 
in  subsequent  sections  we  utilize  Listing  9.2  from  [2]  to  generate  the  von-Karman  phase  screens. 

2.2  Zernike-based  Phase  Screen  Generation 

Drawing  a  stark  contrast  to  Eourier  methods,  the  Zernike  polynomial  approach  does  not  begin  with  a  random  phase 
array,  instead  the  coefficients  are  combined  into  two-dimensional  random  functions.  In  his  1976  paper,  Noll  [6] 
defined  a  new  modified  set  of  Zernike  polynomials  Zy  (p,  6).  He  defined  these  polynomials  to  be 


Zeren  /  =  'Jn  +  cos(m0)] 

, _  ^  >,  m  ^0 

=  Vrmfi™V2  sin(m0)  j  (3) 

Zj  =  Vn  +  m  =  0 


where  the  functions  i?™  (r)  are  referred  to  as  radial  functions  defined  to  be 

(n-m)/2 

nm.  ^  _  V  (-l)"(n-s)! 

Z,  s![(n-tm)/2-s]![(n-m)/2-s]!^  ' 

s=o 

Based  on  a  polar  coordinate  system,  R  is  the  radius,  p  is  the  position  along  the  radius,  and  0  is  the  angle  with  respect 
to  the  x-axis.  NolFs  numbering  sequence  of  the  index  j  corresponds  to  each  polynomial  Zj  and  proceeds  by  row, 
then  given  value  of  radial  degree  n,  and  increases  with  azimuthal  frequency  m.  We  generate  random  phase  screens, 
0atm’  from  the  relationship; 


0atmiRp.S)  =  ^ajZj(p,0). 

j 

To  generate  the  phase  screen  with  accurate  correlation  we  calculate  the  amplitude  ay  vector  by  multiplying  $  with  a 
zero-mean,  unit-variance  random  vector. 


fly  =  (6) 

^  is  the  Cholesky  decomposition  of  the  covariance  matrix 


<P  = 


and 


=  Cy  i 


(7) 


where  the  covariance  matrix  Cjj’  is  generated  from  the  covariance  of  two  Zernike  polynomials  Zy  and  Zy  with 
amplitudes  ay  and  ay; 


r  n  _ K,,S,r[(n  +  n'-5/3)/2](D/rg)^/^ _ 

i.i  J  p|-(„  _  _  i7/3)/2]r[(n  -  n'  -  17/3)/2]r[(n  +  n' -  23/3)/2]  ’ 


(8) 


in  which 


r(14/3)[(24/5)r(6/5)]=/6[r(ll/6)]2 

- ^ 


(_l)(n+n'-2m)/2^(j,  +  l)(n'  +  1)  , 


(9) 


and 


=  (m  =  m)  ^^^parity(J,j')  V  (m  =  0)|. 


(10 


3.  TEMPORAL  CORRELATION 

Without  determining  some  factor  of  correlation  between  individual  pixels,  the  FT  methods  cannot  be  evolved,  only 
shifted  with  statistical  accuracy.  The  benefit  of  the  Zernike  method  however,  is  that  once  a  single  random  screen  is 
generated,  the  correlation  between  each  pixel  is  inherent  within  the  Zernike  polynomials.  Cain's  [4]  work  with 
atmospheric  tilt  has  already  generated  a  statistically  accurate  algorithm  to  temporally  evolve  the  first  two  Zernike 
polynomials.  We  expand  the  algorithm  of  the  first  two  polynomials  to  any  N  maximum  number  of  Zernike 
polynomials,  and  utilize  the  methodology  of  the  expanded  algorithm  to  generate  a  temporal  correlation  of  the 
aforementioned  n  vector.  With  this  temporal  correlation  we  can  evolve  the  a,  vector  with  accurate  correlation  of 
contiguous  phase  screens.  We  utilize  Bayes'  theorem  to  update  the  vector  n  with  mean  Eq.  (11)  and  variance  Eq. 
(12): 


E 


£'[n(t2)|n(ti)]  = 


Var(ji(t^)) 


n(t2)  -n(ti) 


fin(At) 


4 -4m 


(11) 

(12) 


and  solve  for  the  correlation  of  the  n  vector  in  time,  i?^(At).  We  utilize  the  relationship  in  Eq.  (1 1)  to  define  the 
correlation  between  two  it’s  as 


fijj(At)  =  E 


'  N  X  /  N  ^ 


We  then  build  upon  the  definition  of  these  a/s  from  [4]  and  expand  Eq.  (13)  by  substituting  in 


afe(ti)  = 


/ /4(H',s)(p;c(w,s)0Qj,„(w,s)dwds 


a-iih)  = 


/  /  4  (w,  (w,  s)  dwds 

/  /  A(w,s)(pi(w,s)eatmiw  +  (yx)At,s  +  (iVy)At)dw2ds2 
f  f  A(w,s)q)i/(w,  s)dwds 


yielding 


fin(At)  =  E 
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*  ^msi)<Pk(w,s)&atmmsi)dw^ds^ 

r  r  .  y  -.y  vyy.  1^ 
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/  /  A(w,s')(pj^'^(w,s')dwds 

-iff  Mw2,S2)(Pl(W2,S2)0atmiw2+Vx^t,S2+Vy^t)dW2ds- 

f  f  A(w,s)(pi'^(w,  s)dwds 


(13) 


(14) 


(15) 


We  have  redefined  our  notation  of  the  Zernike  polynomials  Zy  as  and  (p;  allowing  for  temporal  correlation. 
Additionally  we  represent  our  aperture  as  A,  and  the  atmospheric  phase  screen  6atm^  ^^ch  with  the  same  coordinate 


system  (tv,i).  By  rearranging,  recognizing  that  the  expectation  operation  is  linear,  and  grouping  similar  terms  we 
redefine  Eq.  (15) 


/  /// A(Wi,Si)A(W2,S2)yfe(Wi,Si)y;(W2,S2) 

f  f  Aiw,s)(p^^(w,s)dwds  j  j  Aiw,s)(pi^(w,s)dwds 
X  ■fi'[0atm('^l.Si)0a£m('V2  +  VxAt,S2  +  V^At)]dWj_dSj_dW2dS2. 


(16) 


This  is  possible  due  to  the  fact  that  the  first  line  in  Eq.  (15)  is  deterministic  given  the  input  parameters  D,  rg,  and  the 
respective  wind  speeds.  Recognizing  that  the  expectation  is  a  correlation  we  are  able  to  redefine  Eq.  (16)  as  Eq. 
(17)  by  a  change  of  variables,  =  W2  —  &  Ty  =  S2  —  Sj^  and  Sx^  =  SW2,  SXy  =  SS2 


RfliM)  = 


N  N 


11 

i'  =  1  7=1 


////  A(Wi,Si)A(wi  +  T^.Si  +  Ty)(pi,(Wi,Si)(pi(wi  +  T^.S^  +  Ty) 

/  /  A(w,  (w,  s)dwds  J  J  A(w,  s)q)i'^ (w,  s)dwds 


X  Rg  [(t;^  +  (v^  +  Vt,)At,  Ty  +  (v^  +  Vg)At)]dwj^dsj^dT^dTy. 


(17) 


We  simplify  the  first  line  of  the  numerator  as  the  independent  multiplication  of  Zernike  polynomials  with  the 
aperture 


F  =  T(wi,Si)(pfe(wi,Si) 

G  =  A(w^  +  T^,S^  +  Ty)(pi(w^  +  T^.S^  +  Ty) 


(18) 


which  we  then  correlate  using  Eourier  transform  properties  and  define  the  inner  double  integral  operation  as  P^  i. 

00  00 

Pk.li'^X.T^y)  =  j  I  A(Wi,Si)¥Jfe(Wi,Si)A(wi  +T;„Si  +T3,)(P;(wi  +T;„Si  +T3,)  dWidSi  (19) 


This  simplifies  Eq.  (17)  to 


—  VV0  “V  -1  +  Vx^t)]dT^dTy 

4712^.^='^  —  f  f  A(w,s)q)i^'^(w,s)dwds  f  f  A(w,s)q)i'^(w,s)dwds' 


(20) 


Unfortunately  we  do  not  have  a  method  for  computing  the  correlation  of  the  atmospheric  phase  screens  Rg  directly. 
Therefore  we  use  the  relationship  of  the  phase  structure  to  the  phase  correlation  function  [4] 


DniAt)  =  2fis(0)  -  2Rfim- 


(21) 


We  rearrange,  add  zero  (+  Rg  [0,0]  —  Rg  [0,0])  and  substitute  back  into  Eq.  (20)  giving 


Off  = 


/  /Tfci(T;„T3,)2 


Rg[0,0]-Rg[0,0]  +  Rg[T^,Ty] 

Reiix^  +  v^At.Ty  +  v^At)] 


dTrdx,, 


(22) 


An'^  =  f  f  A(w,s)(pi^^(w,s)dwds  f  f  A(w,s)(pi^(w,s)dwds 


By  pairing  the  phase  variances  with  the  appropriate  phase  correlations,  recalling  Eq.  (21),  the  structure  function  can 
be  expressed  in  terms  of  a  phase  structure  function  difference. 


-1  .  -1  /  /  Bfc.i(u.Ty)[£>e(T;,  +  V^At,Ty  +  V^At)  -  Dg(T^,Ty)]dT^dTy 
"  An^  f  j  A(w,s)q)i^'^(w,s)dwds  f  f  A(w,s)q)i'^(w,s)dwds 


(23) 


[4]  provides  the  methodology  to  analytically  calculate  the  phase  structure  difference  for  the  aforementioned  inputs. 
We  can  now  solve  for  £>^,  and  by  the  relationship  of  Eq.  (21)  solve  for  7?^.  With  7?^  calculated  for  a  given  N 


polynomials,  D,  rg,  and  wind  speeds  &  Vy,  we  can  accurately  update  the  random  n  vector  with  Eq.  (11)  and  Eq. 
(12).  This  allows  us  to  temporally  evolve  the  phase  screen  (piRp,  9)  with  the  correct  correlation  quickly,  since  the 
Rfl  vector  only  needs  to  be  generated  once  for  any  given  condition.  While  computationally  it  is  inefficient  to 
regenerate  the  vector  each  time  a  simulation  is  performed,  having  a  database  of  stored  vectors  allows  for  very 
efficient  phase  screen  generation  with  the  accurate  correlation.  In  contrast  to  the  frozen  flow  method  in  which  a 
large  phase  array  must  be  generated  and  then  shifted,  this  Zernike-based  turbulent  flow  method  requires  only  the 
correlation  matrix,  C^',  the  correlation  vector,  and  the  Zernike  polynomials,  Zj  to  be  stored.  The  update  then 
becomes  an  extremely  fast  matrix  operation,  and  does  not  require  an  enormous  amount  of  memory  and  processing 
time  to  store  and  shift  giant  phase  arrays. 


4.  ANALYSIS 

To  validate  this  turbulent  flow  approach  to  temporal  evolution  we  focus  on  two  metrics:  temporal  tilt  correlation  in 
time,  and  the  long  exposure  optical  transfer  function  (OTE).  [4]  provides  the  methodology  to  calculate  the 
theoretical  tilt  correlation  given  input  parameters  D,  rg,  and  wind  speeds  8i  Vy,  and  is  based  on  the  accepted  thin 
lens  approximation  phase  structure  of  the  atmosphere.  We  present  this  as  the  theoretical  basis  for  our  x  &  y  tilt  in 
time  analysis. 


,(Cn)  = 


C  C  s)9atm(w,  S,  t„)(p^yCw,  S) 


(24) 


The  tilt  a^y  is  calculated  from  the  first  two  Zernike’ s  cp^y  which  correspond  to  x  and  y  tilt  respectively,  the  aperture 
A,  and  instant  draws  of  our  successive  Zernike  and  FT-based  phase  screens  Since  tilt  is  zero  mean  [4],  to 

calculate  the  correlation  of  these  tilts  in  time  we  recognize  that 

Raitl.W  =  E[ax,yitl),ax,yit2)]  (25) 


and  average  thirty  simulations  to  produce  our  mean  p  measured  tilt  correlations  in  time.  We  performed  this  analysis 
on  four  different  low  wind-velocity  conditions  and  averaged  our  results  to  calculate  the  mean  p  and  standard 
deviation  a^.  We  present  our  simulated  measurements  for  comparison  to  the  theoretical  (underlined)  values  in  Table 

1. 


Table  1.  Monte  Carlo  tilt  correlation  results 


X  Correlation 

Y  Correlation 

von- 
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Zernike 

2 

Theory 
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Zernike 
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Theory 

Karman  p 

A' 

A' 

Karman  p 

A' 

A' 

K  =0.125 

17.873 

5.752 

24.120 

5.862 

26.028 

16.001 

5.589 

25.791 

9.434 

26.0289 

T,  =0.125 

5.971 

5.971 

5.971 

5.971 

K  =0.25 

14.1568 

3.486 

24.653 

6.464 

25.5015 

14.706 

6.088 

26.949 

9.948 

26.0308 

11 

o 

o 

4.655 

4.655 

4.752 

4.752 

V,  =0.25 

13.51712 

4.935 

25.113 

4.745 

24.6239 

10.868 

3.289 

23.123 

2.243 

23.299 

II 

p 

La 

4.495 

4.495 

4.253 

4.253 

K  =0.75 

8.61027 

3.519 

20.454 

2.046 

20.6438 

9.369 

4.071 

22.620 

3.32 

22.2606 

in 

o 

II 

3.734 

3.734 

4.129 

4.129 

This  analysis  clearly  shows  that  the  turbulent  flow  Zernike-based  approach  significantly  outperforms  the  frozen  flow 
von-  Karman-based  approach,  in  regards  to  temporal  tilt  correlation.  The  von-Karman  measurements  provided  for 
each  wind  based  scenario  depict  that  in  all  cases  the  measured  mean  is  off  from  the  theoretical  mean  by  a  factor 


greater  than  one  standard  deviation.  Additionally,  all  Zernike  measurements  fit  well  to  the  theoretical  mean  and  are 
well  within  a  standard  deviation  of  the  theoretical  mean. 

Next  we  examine  the  long  exposure  OTF.  We  generate  the  theoretical  OTF  as  a  product  of  the  diffraction  limited 
OTF  and  Eq.  (26),  [5].  The  diffraction  limited  OTF  is  generated  by  the  FT  of  the  point  spread  function  (PSF) 
created  by  the  pupil  function  with  zero  phase. 
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(26) 


To  simulate  the  measured  long  exposure  OTF  we  simulate  more  than  five  hundred  iterations  of  the  temporal 
evolution,  both  frozen  flow  and  turbulent  flow,  capture  the  instantaneous  OTF,  and  average  the  instantaneous  OTFs. 
The  average  OTFs  are  then  compared  to  the  theoretical  long-exposure  OTF  [5]  and  quantified  via  mean  squared 
error  (MSE)  analysis.  Our  Monte-Carlo  MSE  results  are  provided  in  Table  2. 


Measured  &  Theoretical  OTFs 


Measured  &  Theoretical  OTFs 


(a)  (b) 

Fig.  1:  Long-exposure  OTF  compared  to  measured:  (a)  modified  von-Karman  frozen  flow  OTF,  (b)  Zemike  turbulent  flow  OTF 


Table  2.  Mean-Squared  Error  results 


MSE 

von- 

Karman 

Zernike 

K  =0.125 

Vy  =0.125 

1.17E-04 

2.94E-04 

K  =0.25 

Vy  =  0.0 

2.97E-04 

1.46E-04 

Vy  =0.25 

Vy  =  0.5 

1.18E-04 

9.08E-05 

K  =0.75 
y,  =  0.5 

8.96E-05 

7.51E-05 

In  addition  to  our  findings  on  temporal  tilt  correlation  in  time,  we  see  from  the  long  exposure  OTE  that  this  new 
turbulent  flow  model  is  comparably  accurate  to  the  standard  frozen  flow  model.  Einally  we  note  that  computational 
time  for  each  simulation,  if  treated  independently,  is  comparably  the  same.  Flowever  as  previously  alluded  to,  the 
Zernike-based  approach  described  in  this  paper  allows  for  to  be  calculated  once  for  any  given  conditions.  With 


the  Rfi  vector  already  populated,  the  turbulent  flow  model  produces  temporal  phase  screens  considerably  faster  than 
frozen  flow  model  does. 


5.  CONCLUSION 

Through  the  analysis  provided  above,  we  have  generated  a  method  of  turbulent  flow  phase  screen  generation.  We 
have  shown  that  a  random  phase  screen  generated  in  a  conventional  manner  can  be  temporally  evolved  for  low  wind 
velocity  situations  using  this  Zernike  polynomial-based  approach.  This  temporal  evolution  of  phase  as  the  screen 
slowly  shifts  has  been  validated  through  statistical  analysis.  We  compared  this  approach  to  a  widely  used  frozen 
flow  approach  and  highlighted  the  turbulent  flow's  significant  outperformance.  Finally,  having  developed,  tested, 
and  validated  this  turbulent  flow  model,  we  can  perform  further  analysis  on  the  atmospheric  effects  of  ultra  long 
pulse  LADAR  systems. 
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